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Abstract. In many clusters of galaxies there is evidence for cooling flows in the central 
regions. The ultimate fate of the gas which cools is still unknown. A possibility is that a 
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to H2, HD and CO molecules. As an example, we determine the minimum temperature 
, achievable by clouds in the cooling flows of the Centaurus, Hydra and PKS 0745-191 
. clusters. Our results suggest that clouds can reach very low temperatures - less than 

~ 10 K - which would explain the non-detection of CO in these clusters. 

Key words: galaxies: clusters - galaxies: cooling flows - X-rays: galaxies - molecular 
processes 



1. Introduction 

In several clusters of galaxies one observes a soft X-ray excess towards the central 
regions. This excess is interpreted as being due to hot intracluster gas (10^-10^ K) 
with a cooling time less than the Hubble time (Fabian 1994). This gas cools and falls 
quasi-hydrostatically into the center of the cluster potential well (e.g. Fabian, Nulsen & 
Canizares 1984, 1991; Sarazin 1988). 
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Indeed, radial profiles of temperature and density in the intracluster gas inferred from 
X-ray observations show that cooler denser gas (10^-10^ K) is present in the central 
region of the cluster (Johnstone, Fabian. Edge & Thomas 1992). Typical mass accretion 
rates inferred in these cooling flow clusters are ^ 100 Mq year~^ and even higher in some 
cases (Fabian 1994). Thus cooling flows in clusters of galaxies deposit large quantities 
of cool gas around the central galaxy, which is still continuing to grow today. Recent 
extreme ultra-violet (EUV) observations revealed the presence of important amounts of 
a warm EUV emitting gas in the cooling flow cluster A1795 (Lieu, Mittaz, Bowyer ct al. 
1996a, 1996b). 

The final evolution of the cool gas is not clear. It may just accumulate as cool dense 
clouds, indeed White. Fabian, Johnstone ct al. (1991) have observed a large amount of 
cold X-ray absorbing matter distributed over the inner few 100 kpc of some clusters of 
galaxies and came to the conclusion that the absorbing component is likely to be either 
atomic or molecular (mainly hydrogen). Moreover, they suggested that the gas is very 
cold (< 10 K) and that it makes up a significant fraction of the total mass, which has 
cooled out of the cooling accretion flow during the lifetime of the cluster. 
Fabian, Johnstone & Daines (1994) have postulated that dust can form in the cold clouds 
embedded in cooling flows, which would explain why the clouds are almost undetectable 
outside the X-ray wavelength range. However, the nature of the absorbing material re- 
mains uncertain. The limits imposed by 21-cm and CO observations, suggest that it is 
unlikely that signiflcant amounts of cold gas can remain undetectable. On the other hand 
Jaffe & Bremer (1997), using K-band spectroscopy, have found that the inner few kpc 
of central cluster galaxies in cooling flows have strong emission in the iJ2(l-0)S(l) line. 
This is not seen in a comparable sample of non cooling flow galaxies. Jaffe & Bremer 
(1997) suggest that it is likely that a large mass of additional molecular material is in 
the cooling flow, producing the soft X-ray absorption. At least a fraction of the cool gas 
must reach the inner regions of the central cluster galaxies, where it could accumulate 
into molecular clouds and subsequently form stars (Cardiel, Gorgas & Aragon-Salamcnca 
1998). 

The metallicity of a cluster seems to be correlated with the presence of a cooling flow 
(Allen & Fabian 1998). Mushotzky & Lowenstein (1997) showed that most of the en- 
richment of the intercluster medium occured at high redshifts. This is consistent with 
the semi-analytic models of galaxy formation (Kauffmann & Chariot 1998), which find 
indeed that the metal enrichment occurs at z > 1. In this context molecules such as CO 
can subsequently be formed in the gaseous medium. The molecular gas can cool down 
to very low temperatures within the intra-cluster environment as has been suggested by 
Ferland, Fabian & Johnstone (1994). 
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O'Dea, Baum, Maloney et al. (1994) searched for molecular gas, by looking for CO emis- 
sion lines in a heterogeneous sample of five radio-loud galaxies (three of which are in 
cooling flow clusters), using the Swedish-ESO Submillimeter Telescope in the frequency 
ranges of 80-116 GHz and 220-260 GHz. A positive detection of CO has been made in 
two of the galaxies: PKS 0634-206, a classical double radio galaxy, and possibly in NGC 
4696. They come to the conclusion that in order to have escaped detection the gas has 
to be very cold, close to the temperature of the Cosmic Background Radiation (GBR). 
O'Dea, Baum, Maloney et al. (1994) estimated also the expected temperature of the 
clouds by comparing the X-ray heating with the cooling rate for molecular gas. For the 
latter they adopted the cooling rate calculated by Goldsmith & Langer (1978), which 
includes line cooling from both molecular and atomic species likely to be abundant in 
molecular clouds. They thus assumed that the clouds have a metallicity similar to solar 
neighbourhood clouds. However, the Goldsmith & Langer cooling rates are valid in the 
range ~ 20 - 60 K, below ^^ 20 K they are no longer reliable. In addition the possible 
role of HD as a cooling agent was not taken into account. Therefore, the conclusion of 
O'Dea, Baum, Maloney et al. (1994), that the minimmir temperature of the clouds is in 
the range 20-30 K is certainly questionable. 

Recently, Bridges & Irwin (1998) have observed molecular gas in the Perseus cooling flow 
galaxy NGG 1275 and have mapped the central arcminute of this galaxy with the James 
Glerk Maxwell Telescope in ^^C0(2 — 1). Assuming a galactic conversion between CO 
intensity and the molecular H2 gas mass they inferred a total mass M/f^ ^ 1.6 x 10^" Mq. 
They argue that the reason why CO has been detected in NGG 1275 is that the radio 
emission heats up the molecules to a detectable temperature. Indeed, NGG 1275 happens 
to have one of the strongest radio core of any cooling flow central galaxy. Without the 
heating of the radio emission the gas would stay at much lower temperatures and thus 
escape detection. 

The aim of this paper is to investigate the minimum temperature achievable by clouds in 
cooling flows by improving the analysis done by O'Dea, Baum, Maloney et al. (1994), in 
particular by considering clouds made of H2 , HD and CO molecules and by computing 
cooling functions which arc more appropriate for temperatures below 20 K. We do not 
discuss the different pathways for forming H2, HD and CO molecules. H2 and HD can 
be formed at the early epoch of the recombination of the hydrogen (see Lepp & ShuU 
1984, Puy, Alecian, Lc Bourlot et al. 1993), whereas CO formation can be triggered by 
the radiative association of C with O after the first star formation. 
The paper is organised as follows: in Sect. 2 we recall the effects of X-ray heating on the 
clouds using the results of Glassgold & Langer (1973) and O'Dea, Baum, Maloney et al. 
(1994). In Sect. 3 we calculate the molecular cooling due to iJ2, HD and CO, respectively. 
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through the transition between the ground state and the first rotational level, since we 
are interested on what happens for very cold clouds (with temperature of order 20 K or 
even below). Moreover, also from the observational point of view the transition between 
the ground state and the first excited state gives possibly the strongest contribution at 
very low temperatures as suggested by the measurements by Jaffe & Bremer (1997) on 
the emission in the H2{l-0) S(l) line. We then estimate in Sect. 4 the equilibrium tem- 
peratures of the clouds for the Ccntaurus, Hydra and PKS 0745-191 clusters, and discuss 
the observational consequences. 
Some technical details of the calculations are given in two Appendices. 

2. X-ray heating of the clouds 

The clouds are embedded in the hot intracluster gas, whose emission is dominated by 
thermal bremsstrahlung and thus the clouds are mainly heated by X-rays. The mecha- 
nisms by which soft X-rays heat a gas cloud have been extensively studied by a number 
of authors: Habing & Goldsmith (1971) and Bergeron & Souffrin (1971) for Hj regions, 
Glassgold & Langer (1973) for molecular hydrogen clouds. 

Due to the increased number of possible charged-particle reactions with molecules, the 
heating mechanism for molecules is different and more complicated as compared to the 
heating processes for atoms alone. Glassgold & Langer (1973) have given a survey of the 
physical processes for a pure molecular hydrogen medium. While X-ray heating of molec- 
ular hydrogen regions has been extensively studied, this has not been the case for other 
molecular species. Glassgold & Langer (1973) discuss in detail the energy loss processes 
and present a quantitative analysis which shows which part of the incident X-ray energy 
will go into heat. In Appendix A we recall the heating processes for H2, HD and CO 
molecules which are of interest for us. 

In order to evaluate the bremsstrahlung fiux coming from the intracluster gas, we must 
define a density distribution for the electrons in the gas. A classical profile is a non- 
singular isothermal sphere (Jones & Forman 1984): 

7ie(i?) = ° ^ 2 cm"^, 

' + (*) 

where no ^ 10^^ cm^'^ is the central gas number density, Vc is the core radius and R is 

the distance from the center of the cluster. In cooling fiow clusters the central density is 

somewhat higher in the range lO^'^ — 10^^ cm^'^ (see Stewart, Fabian, Jones & Forman 

1984). 

For temperatures T > 2 keV the emission from the hot gas is dominated by thermal 

bremsstrahlung and the emissivity for an ion with charge Z and number density Ui in a 
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0.01-0.5 


1.2 X 10-28 


0.1 


-0.65 


0.5-2.5 


4.6 X 10"^^ 


1/2 


-1.1 



Table 1. Parameters for the X-ray heating (sec O'Dca, Baum, Maloney et al. 1994). 



gas with electron density n^ is 



2.4 X 10-'"^ Z^n,n,T^J^'exp{-Ekev/Tkev)9ff 
e°Z'^?ien, (inkeVcm"'"' s~^ kcV"^ sr'^), 



(1) 



where g^^ is the velocity-averaged Gaunt factor, while gas temperature T and photon en- 
ergy E arc given in kcV. The Gaunt factor for the relevant parameter range is accurately 
given by 

TkeV 



gjf - 0.551n(2.225) 



Ek 



eV 



(2) 



The bremsstrahlung "flux" (i.e. 4:TtJe, where Je is the mean intensity at energy E), seen 
by a cloud located at a radial distance R from the center of the cluster is given by O'Dea, 
Baum, Maloney ct al. (1994) in keV cm^^ s^^ kcV"^ 



^irJEiR) 



l.bnnlr'^ 



^^ r r-2fiiTiflrr-2 -^ fi-2^3/2 

dr dO - 
Jo (r2+r2)2(r2 + i?2_2rii!cos0) 



(3) 



{rl + i?2)3/2 

The factor 1.5 takes into account the presence of He and heavy elements. Rb is the outer 
boundary of the cluster gas distribution. The integral in Eq.(3) for typical values r^ = 100 
kpc and Rb ~ 500 kpc with < i? < 500 kpc varies between 1.57 and 7. A reasonable 
average value for the integral is tt, which for clouds located in the inner regions of the 
cluster is an upper bound. 

The problem is much more complicated if wc consider the clastic scattering of electrons 
from H2 molecules. We have various dissociation and ionization channels which contribute 
to the heating (sec Glassgold & Langer 1973). O'Dea, Baum, Maloney et al. (1994) assume 
that the fraction of primary photoclectron's energy which goes into heating is close to 
half. 

Moreover, O'Dea, Baum, Maloney et al. (1994) consider a possible attenuating column 
density within a cloud, which they characterize by the parameter t. Finally, with these 
different approximations they obtain for the heating rate (in erg cm'^ s"^): 



Tx{R) = lo 



nH2n_3 rioo 

572 
1 ■ -" 



Ti 



keV 



(4) 



where the parameters are given in Table 1 and ri_3 = rio/lO ^ cm ^ and rioo — Tc/IOO 
kpc. 
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3. Molecular cooling 

In the standard Big Bang models primordial chemistry took place around the epoch of 
recombination. At this stage the chemical species essentially were H , iJ+, D, D+, He 
and Li. Then with the adiabatic cooling of the Universe due to the expansion, different 
routes led to molecular formation (Lepp & ShuU 1984. Black 1988, Puy, Alecian, Le 
Bourlot et al. 1993, Galh & Palla 1998). 

Heavier molecules than H2 play an important role in the cooling of the denser clouds 
(Lepp & ShuU 1984, Puy & Signore 1996, 1997, 1998a, 1998b) via the excitation of 
rotational levels. Ferland, Fabian & Johnstone (1994) have investigated the physical 
conditions within molecular clouds, particularly the structure of gas clouds in cooling 
flow and showed that molecular cooling plays an important role. 

We calculate the molecular cooling in the range of temperatures ^ 3-120 K, where 120 K 
corresponds to the temperature of the transition between the ground state (J = 0) and 
the first rotational level (J = 1) for the HD molecule. In this range, molecules such as 
i?2, HD and CO are very efficient for cooling. HD and H2 arc the main cooling agents 
between 80-120 K, then HD dominates between ~ 40-80 K {H2 cooling, in this case, 
turns out to be negligible). CO can be the main coolant in the range of temperature 
^ 3 K up to '^ 40 K depending on its abundance. In order to calculate analytically 
the molecular cooling, we only consider the transition between the ground state and the 
first rotational level, although below 120 K it is well possible to excite other rotational 
levels of CO, for instance the next higher transition (J = 1 ^ J = 2) corresponds to 
a temperature of ~ 11 K. The reason is that we are interested on what happens with 
very cold clouds (with temperature of order 20 K or less) and by considering the higher 
levels the molecular cooling gets even more efficient. Thus by taking into account only 
the first level we will get a lower bound on cooling and, therefore, an upper bound on 
the equilibrium temperature. 

The cooling functions for the molecules iJ2, HD and CO arc computed in Appendix B. 
We have taken into account radiative transfer effects following the treatment by Castor 
(1970), Goldreich & Kwan (1974) in a similar way as done in the paper of Goldsmith & 
Langer (1978). Thus, the rates for absorption and stimulated emission are given entirely in 
terms of two locally determined quantities, the source function and the escape probability 
Po which depend on the molecule (see Appendix B). Finally we obtain for the cooling 
functions for H2, HD and CO (in erg cm^^ s~^): 

Ah, - 8.61 X 10-^4 P^^n^j^ VT™ ^^P^-f^l^ra) ^ ^j.) 



' 0.69/3?^ +V7T^nff, 



1 + 5ea;p(-512/r„ 



A 0Rr^in-2i«i/n 2 f7^ ea;p(-128.6/T„j) 

KhD ^ 2.66x10 13^ VHOlT-H^y/Tm ^^ ^ f ■ ■ — ,[&) 



Uiep^D + Vr^riH^ l + 3exp{~128.6/T, 
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with the ratios tjhd = nHD/nH2 a-nd "^co ~ "■co/"--ff2 • The various quantities entering 
in the above equations are defined in Appendix B. 

Combes & Wiklind (1998) have detected molecular absorption lines systems with the 
15-m Swcdish-ESO Submillimetre Telescope at high redshift. In particular, they have 
detected CO absorption lines in some galaxies or quasars. For Centaurus A they estimated 
the column density for CO and H2 to be Nco = 10^^ cm~^ and N/j^ ~ 2.0 x 10^" cm~^, 
respectively. For some clusters they could only derive an upper bound of about 6 x 10^^ 
cm"^ for the CO column density. The values for the CO and H2 column densities can 
vary substantially between different clusters. 

Moreover, we do not know if the column density is due to the sum of several smaller clouds 
which are present along the line of sight or to a unique big cloud. It is thought that a 
Jean's unstable and isothermal cloud fragments into a small number of subclouds, that 
in turn repeat the process recursively at smaller scale as long as isothermal conditions 
prevail. Thus it is not unrealistic to consider very small clouds in cooling flow as a result 
of fragmentation (Fabian 1992). Indeed, the high pressure in a flow implies that the Jeans 
mass of a very cold cloud is small (Fabian 1992). We expect than that the molecules are 
mostly produced in the small clouds where the densities can be higher with respect to 
the surrounding regions and therefore also the cooling will be most efficient there. 
These small clouds are most likely embedded in bigger less dense clouds with decreasing 
molecular abundance and also higher temperatures. However, the gas surrounding the 
small clouds will attenuate the incoming bremsstrahlung flux as taken into account by 
the T parameter discussed in Section 2. 

A scenario with successive fragmentation of the clouds in cooling flows has been proposed 
by Fabian & Nulscn (1994) and in another context, namely for the dark matter in the 
disk of a galaxy by Pfenniger & Combes (1994), whereas in the halo of our galaxy by De 
Paolis et al (1995). Pfenniger & Combes (1994) have proposed that the cold H2 clouds 
have a fractal structure. 

We will thus in the following assume as a picture that the clouds fragment into small 
ones as a result of cooling processes. For definiteness we will adopt the following column 
densities for a typical small cloud: Nco ~ 10^"* cm^^, iV^a = 2 x 10^^ cm^^ with 
nH2 = 10^ cm^'^ (density of H2 the collisional species). We have the following ratio 
Vco ^ j^ ~ 5 X 10^^. O'Dea, Baum, Maloncy et al. 1994 instead use the value 
r]co ~ 10^^. For HD we take r]HD ~ %"" ^ 7 x 10^^, which corresponds to the 
primordial ratio (Puy, Alecian, Lc Bourlot et al. 1993). This should possibly be an upper 
limit. These characteristics correspond to a size for the cloud of L ~ 10^^ pc which is 
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T„ 


2.735 


5 


20 


10 


/?f° 


0.47 


0.72 


0.89 


0.96 



Table 2. Escape probability f]^'^ for a photon emitted in a transition J = 1 ^ J = 
for the CO molecule, assuming a column density Nco = lO^"*^ cm^-^ of the cloud. 

close to 30 AU as for instance proposed in the fractal model of Combes & Pfenniger 
(1994). 

Using the above values we can then estimate the escape probability for a photon emitted 
in a transition in the temperature range ^ 3 K to 120 K due to the three considered 
molecules. The escape probabilities due to HD and H2 are close to 1 {p^^ = (3^^ ^ 1) 
in the range where HD and H2 are dominant {T„i > 40 K). In the case of CO molecule, 
which is supposed to be the main cooling agent at very low temperature, the probability, 
that the photons emitted during the transitions J = 1 ^ J = are reabsorbed, is non 
negligible and thus the escape probability (3^'^ is below 1. In Table 2 we give some values 
of P^'-' for different temperatures. Particularly at low temperatures re-absorption has to 
be taken into account and leads thus to a less efficient CO cooling. 

Our calculations differ from the ones obtained by O'Dea, Baum, Maloney et al. (1994), 
since they used the cooling rates calculated by Goldsmith & Langer (1978), which in- 
clude line cooling from both molecular and atomic species likely to be abundant in disk 
molecular clouds. Moreover, the possible role of HD as a cooling agent was not taken 
into account. Goldsmith & Langer (1978) considered molecular hydrogen densities in 
the range 10^ to 10^ cm^'^, kinetic temperatures between 10 K and 60 K and a wide 
variety of molecules with fractional abundances bigger than or comparable to 10~^ (^2, 
CO, H2O, O2, C2, N2 and hydrides). They find, hke us, that below 40 K CO is the 
dominant coolant. Nevertheless, they approximated the cooling rate by a power law for 
temperatures between 10 K and 60 K. Instead, we approximate the molecular cooling by 
considering only the transition between the ground state and the first rotational level. 
We examine the equilibrium of clouds for temperatures below 10 K, a range where the 
power law dependence of Goldsmith & Langer (1978) is no longer valid. 

4. Results and discussions 

We have seen that molecular clouds in cooling flows are heated by the external X-ray 
flux. At equilibrium the heating and molecular cooling arc equal: Tx = A^, where the 
total coohng function is defined as At = A^^ + Ahd + ^co, see Eqs. (5), (6) and (7). 
The A-ray heating rate (see Eq. 4) depends on the hot gas parameters n_3, rioo and 
temperature T^ev, for which we adopt the values (see Table 3) of three clusters considered 
by O'Dea, Baum, Maloney et al. (1994): Hydra A, Centaurus and PKS 0745-191. 
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Cluster 


n_3 


rioo 


TkeV 


References 


'^cool 


References 


Hydra A 


6.5 


1.45 


4.5 


David et al. 1990 


162tll 


Cardiel et al. 1998 


Centaurus 


9 


2 


2.1 


Matilsky et al. 1985 


871L 


Allen & Fabian 1997 


PKS 0745-191 


35 


0.5 


8.6 


Arnaud et al. 1987 


214I25 


Cardiel et al. 1998 



Table 3. Parameters for the galaxy clusters we consider in the text, n_3 is in units of 
10~^ cm^'^, rioo = rc/100 where re is in kpc and TkeV in keV (see O'Dca, Baum, Maloney 
et al. 1994). rcooi is the cooling radius with upper and lower limits. 



Tl 


T2 


T-3 


r4 


Ts 


T"6 


0.01 


0.25 


0.5 


1 


1.5 


2.5 



Table 4. Values of the attenuation factor t, which will be used in the text. 



We have seen (Eq. 4) that Tx depends on the distance R (from the center of the 
cluster) at which the molecular cloud is located. Thus at the equilibrium temperature 
Teg, where T{Req) ~ KT{T(.q), we obtain 



^eq ^c 



-'jon'i^nH^rwoTl^ev'^^ 



Al 



2/3 



-1 



(8) 



The intracluster gas is denser in the core of the cluster and, therefore, the radiative 
cooling time, tcooh due to the emission of X-rays is shortest there. A cooling flow is 
formed when tcooi is less than the age of the system t^ysi roughly tgyg ^ H~^ where Ho 
is the Hubble constant. We define the radius rcooi where tcooi ^ ^sys, beyond this radius 
we do not expect a cooling flow. Table 3 gives the cooling radius for the three clusters 
we consider. 

In Figs. 1, 2 and 3 we plot the curves Req as a function of T^q for different values of 
the attenuation factor r (sec Table 4) for the AT-ray heating, which is due to the presence 
of an attenuating column density (O'Dca, Baum, Maloney et al. 1994). Upper curves in 
the different figures correspond to low values of t, in which case the AT-ray heating is 
more important. Thus low equilibrium temperatures are achieved only at large distances 
from the cluster center. We indicate on the plots also the cooling radius rcooi- We notice 
that the curves above rcooi in the various figures are given just for clarity, but we do 
not expect them to be realistic, since the density of the clouds will decrease in the outer 
regions of the cluster, so that our adopted values for the density will be too high and thus 
no longer valid. The slope of the Rcq curves depends on the molecular cooling At (see 
Eq. 8). For all cases, H2 cooling turns out to be negligible for equilibrium temperatures 
below 80 K. CO cooling is dominant in the range from ~ 3 K up to 40 K, whereas HD 
dominates between ^ 40 K and 80 K. 
Fig. 1 is for the Hydra A cluster. Thermal equilibrium in the cooling flow region 
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{Req < Tcooi) IS posslble for different values of r. Notice, that for t = T5 and t = tq 
the equilibrium temperature T^q can be as low as ^ 10 K. If R^q ^ 200 kpc and rcooi 
close to the upper limit (218 kpc), it is even possible that Te,q is as low as ^^ 3 K for rg. 
O'Dea, Baum, Maloney et al. (1994) suggested, in order to explain the non-detection of 
CO in this cluster, that the gas is molecular and very cold with temperature close to 
that of the CBR. Our calculation seems to confirm this conclusion. 
Fig. 2 is for the Centaurus cluster, for which the situation is quite different, since the 
cooling radius is smaller (~ 87 kpc). For this reason the equilibrium between molecu- 
lar cooling and X-ray heating is difficult to realize at small radii. We find equilibrium 
solutions in the cooling flow region (not for ti however) for a temperature beyond 40 
K, for which HD cooling is important. In this case the temperature of the clouds is 
not sufficiently low in order to explain the non-detection of CO, unless rjco is higher 
or the total gas mass is relatively small. This latter possibility has also been mentioned 
by O'Dea. Baum, Maloney et al. (1994) based on their observational upper limit for the 
J = ^ ,7 = 1 transition in CO for NGC 4696. 

In Fig. 3 we consider the PKS 0745-191 cluster. Although the cooling radius TcooI — 214 
kpc is the largest one, the temperature T^eV ~ 8.6 keV of the cluster is very high and 
thus the X-ray heating is very important. Therefore, for this reason the equilibrium be- 
tween molecular cooling and X-ray heating for low temperatures would be possible only 
at large distances. In the region of cooling flow {Req < Tcooi), we find that the equilibrium 
temperature is possible at very low temperature {T^q ~ 5 K) for T5, Tg and at Tf,q ^ 3 
K for Tg only. Thus, as for the Hydra A cluster, this could explain the non-detection of 
CO. 

Of course it is difficult to draw firm conclusions on the existence of very cold clouds in 
cooling flows, since an important parameter like the CO abundance is quite uncertain. In 
Figs 4, 5 we plotted the curves Req (for Hydra A and PKS 0745-191) as a function of Teq 
for different values of r/co (see Table 5) keeping n/fj = 10^ cm^-^, Nh2 = 2 x 10^^ cm~^ 
and T = tq fixed (keeping thus the same size L ^ 10^^ pc of the cloud and changing 
Nco)- From the figures we see that at very low temperatures Req varies by about a 
factor of 2 when increasing ^)co 10 times. If we took into account all transitions, in 
particular for r,„ > 10 K, the dependence on the CO abundance would turn out to be 
less important. We notice that there are other values of rjco for which it is possible to 
have an equilibrium temperature close to the CBR one. If Nco is substantially smaller 
than 10^"* cm^^ ('7co << ilco,i) the clouds are optically thin but cooling gets less 
efficient. The net effect will be an increase for the value of Req- On the other hand, if 
Nco is bigger (i.e. r]co > Vco.s) we have only a small decrease of Req at Teq ~ 3 K. 
The hypothesis of very cold molecular gas in cooling flows seems reasonable given also 



Puy, Grenacher, Jetzer: Thermal Equilibrium.. 



11 



300 



200 






K 



-n 



100 



I'M 

\ \ 

\ \ 
\ 



HYDRA A 



\ 



W 



\ \ 



\ \ \ \ 



h 1 \- 

\ 

-V 



X \ \ \ 



\ 



\ \ \ \ \ 

\ \ 



T>. 



20 



40 



60 
Tea [K] 



80 



100 



120 



Fig. 1. Radius R^q in kpc as a function of the cloud temperature T^q for the Hydra 
cluster. Upper curves correspond to low r values. The flat curves indicate the cooling 
radius (mean value with upper and lower limit), rjED = 7 x 10^^, rjco = 5 x 10^^ and 
nH2 ~ 10^ cni^"^. For ti, T2 and T3 we take for 70, a and S the values given in the first 
line of Table 1. whereas for t^, t^ and rg the corresponding values are given in the second 
line. 



the fact that with our approximations we get upper limits for the cloud temperatures. 
We notice, that two clusters have an important cooling flow: 400 Mq yr~^ for the Hydra 
cluster (David, Arnaud, Forman & Jones 1990) and 500 Mq yr^^ for the PKS 0745 
cluster (Arnaud, Johnstone, Fabian et al. 1987), whereas the Centaurus cluster has a 
smaller cooling flow of 20 M© yr^^ (Matilsky, Jones & 1985). An important cooling flow 
could trigger thermal instabilities in the gas (David, Bregman & Scab 1988, White & 
Sarazin 1987) and could lead to the fast-growing of perturbations, which can then form 
objects such as molecular clouds. 

Other molecules than CO could also be important, for example CN or H2CO molecules. 
For CN the corresponding temperature for the transition ,7 = -^ J=lis~ 5.44 K 
and for H2CO the transition is close to 6.84 K (see Partridge 1995), so both molecules 
could play a role in cooling. The study of the chemistry in cooling flows could lead to 
important insight in the phenomenon and be of relevance in order to correctly estimate 
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Table 5. Values of the ratio r]c'o, used in the text. 



the total involved gas mass. Such a detailed study, however, is beyond the scope of the 
present paper. 



Puy, Grenacher, Jetzer: Thermal Equilibrium.. 



13 



300 



200 



Hl\ 






K 



100 






PKS 0745-191 



..\. 



il \ ^, '"cool 

H — "^-^ 



v._ \...\. 



\. 



\ 



\ 



\ \ \ \ ^. 



-3 --., 



\ 



A 



J \ W—I I L_ 



20 



40 



60 
Te, [K] 



80 



100 120 



Fig. 3. Radius Req in kpc as a function of the cloud temperature T^q for the PKS 0745- 



191 cluster. See also the caption in Fig. 1. 



14 



Puy, Grenacher, Jetzer: Thermal Equilibrium., 



o 
ft 



K 



ouu 


\ 
V "\^ 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 
HYDRA A 




200 


\\\ 


cool 


100 


" \\\ 

- W "x 


^^~--^^_^ ''^C0,2 


— ' 


n 


''7co,5 

1 1 1 1 


\ 1 1 '"^ 1 
1 1^ 1 1 1 1 1 1 1 t 1 1 1 1 1 





10 12 14 



Te„ [K] 



Fig. 4. Radius Reg in kpc for the Hydra cluster for different values of 7]co {t 
Nh2 = 2 X IQis cni-2 and hh^ = 10^ cm-^). 



T6, 



Puy, Grenacher, Jetzer: Thermal Equilibrium., 



15 



300 



200 



o 
ft 



K 



100 



- \^ 
\ 


1 1 

PKS 


1 1 
0745- 


1 1 1 
-191 


1 1 






- 


\ 
- \ -""cool 


~- — - 


-- — - 


— 






^CO.l 




1^ \ 














- 




1 1 


1 1 


1 1 1 


1 1 


Vc 


^C0,2 
0,3 


- 




- 



10 



12 



14 



Te„ [K] 



Fig. 5. Radius Req in kpc for the PKS 0745-191 cluster for different values oirjco {t ~ tq, 
Nh^ = 2 X 10^* cm~2 and hh^ = 10^ cm-^). 



Acknowledgments 

We thank Monique Signore and Markus Strassle for valuable discussions. We are also 
grateful to the referee for his comments. This work has been supported by the D'' Tomalla 
Foundation and by the Swiss National Science Foundation. 



16 Puy, Grenacher, Jetzer: Thermal Equilibrium... 

Appendix A: Molecular heating 

In this appendix we briefly recall the heating processes for H2, HD and CO molecules, 
which occur mainly through electronic excitations. 

H2 molecules 

The photoionization process of H2 creates an energetic electron and an H2 ion (Wishart 
1979): 

Xray + H2 ^^ H+ + e- . (A.l) 

Different routes are then possible: one of the ways is that iJj recombines with a thermal 
electron (Nakashima, Takanagi & Nakamura 1987): 

H+ + e- ^^ H2 + hiy, (A.2) 

producing a photon but no heating, since the photon can escape from the cloud . The 
second route is that iJ^ collides with H2 and gives rise to the reaction (Prasad & Huntress 
1980) 

H+ + H2 '-^ H+ + H , (A.3) 

followed by the reaction 

H+ + e- ^^ H2+H . (A.4) 

These two reactions are exothermal and lead, therefore, to a transfer of energy and a 
heating of 10.96 eV for every H2 involved in the process (Glassgold & Langer 1973). 

HD molecules 

The fact that the HD molecule has a permanent dipole moment {hhd = 8.3 x 10~^ 
Debye, see Abgrall, Roueff & Viala 1982) leads to X-ray heating processes, which are 
different from those for H2- Thus, after the classical photoionization process (Von Bush 
& Dunn 1972): 

Xray + HD < — > H D+ + e", (A.5) 

one could expect a chemistry similar to reactions (A3) and (A4) with an intermediate 
species H2D~^, that takes over the role of the H^ ion. The knowledge of the i?2-D+ 
chemistry is, however, very poor and we will not discuss it in more details. The most 
exothermical reaction is 

HD+ + e" I — > H + D, (A.6) 

which gives 10.93 eV for every HD^ (see Nakashima, Takanagi & Nakamura 1987). 
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CO molecules 

CO has a permanent dipole moment {^co = 0.118 Debye, see Lang 1974). Exothermal 
chemistry relies on the ion HCO^ (Prasad & Huntress 1980, Anicich & Huntress 1986). 
After photoionisation which produces CO"*": 

Xray +C0 ^^ C0+ + e", (A.7) 

further collision with H2 molecules can produce the intermediate species HCO^: 

CO+ +H2 ' — > HCO+ + H. (A.8) 

Afterwards the dissociative recombination 

HCO+ + e- I — ^ CO + H (A.9) 

gives 7 eV for each process, which goes into heating. 

Electrons 

For the electrons the situation is different. The energetic primary electrons lose their 
energy by electronic excitation and ionisation, the latter process leading to secondary 
electrons: 

e- +H2 ' — > e^ +62 + H+ . (A.IO) 

Whenever the kinetic energy of a primary or a secondary-generation electron falls below 
the threshold for inelastic scattering, it can lose energy only by elastic scattering on 
thermal electrons and thus most of his energy goes into heating. 



18 Puy, Grenacher, Jetzer: Thermal Equilibrium... 

Appendix B: Molecular cooling 

In this Appendix wc compute the coohng functions for the various molecules. 



HD molecules 

The excitation of the ground state is radiative or collisional. The equations of evolution 
for the populations {X^^ for the ground state J = 0, X^^ for the first excited state 
J — \) are given by: 

dX^D 



dt 

dxr 

dt 
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(B.l) 
(B.2) 



where C^^* and C^i^ are the collisional coefficients, which characterize the collisions 

between HD and H2. In the approximation of a Maxwell-Boltzmann distribution we 

have: 

-'i.o 






with T/^q° = Tf^^ = 128.6 K being the transition temperature between the states ,/ = 
and J = 1, T„i the matter temperature of the cloud, B^^ and -Bf^(P the Einstein 
coefficients. The energy density of the cosmic background radiation (T^ = 2.735 K) at 
the frequency v^^ of the transition J = O^J=lis 



HD 

*i,o 
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HD 
^0,1 
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^1.0 



87r/i(z.f^^3 
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/ rpHD 

expi-^ 



512 TT kg Bjjjjfi^j-, 



ih^c^ 3 

is the Einstein coefficient for the spontaneous emission (sec Kutner 1984) and B^d = 
64.3 K is the rotational constant (see Herzberg 1950). 

We have the following normalisation for the populations: X^^ +X{^^ = 1. We consider 
a quasi-instantaneous transition, therefore, we can set 

dxr dxr 

dt dt 

With these approximations we find for the populations X(^^ and X^^: 
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(B.3) 



(B.4) 
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The probabiUty for radiative de-excitation is 

AHD I TDHD..HD 
^1,0 + ^1,0 "l,0 



pHD ^ 



'-'1,0 +^1,0 + -°i,o "l,0 
The molecular cooling function (coUisional excitation followed by a radiative de- 
excitation) is given by 



I^HD — nnoXg Cgi Pj. E 



1,0 ' 



(B.6) 



where E^^ = ksTl^ff is the energy for the transition ,/ = 
With the following relation between the Einstein coefficients: 



B^f = 3B^f and B?.^ 
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we get for the HD cooling function 
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Since T^^ » Tr we can use the following approximations: 
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In this way we obtain for the HD cooling function: 
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This result is valid for the optically thin limit. In the optically thick case one has to take 
into account radiative transfer effects. Goldreich & Kwan (1974) computed the average 
escape probability for a photon emitted in a radiative transition from the J -I- 1 to the J 
rotational level, denoted by f3j, using Castor's result (Castor 1970): 



Pj 



1 - exp(-rj) /tj 



(B.9) 



where tj is the optical depth for the transition J ^ ,/ -|- 1. In our case we consider only 
the transition between the ground state and the first excited level, and thus we get 
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and the corresponding average escape probability is 
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Nhd is the column density of HD. The linewidth Au is assumed to be given by thermal 
doppler broadening 
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(B.12) 
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We neglect the other possible sources of broadening such as turbulence or large-scale 
systematic motions. 

The j3^^ factor enters only in the spontaneous emission term. Therefore, the cooling 
function in Eq. (B8) gets modified as follows 



A 



3Ci«o^ exp{~ ^ ) hv^^^ uhdA^^HI 



HD 



HD 



A^.^P?"" + Cfo^ 1 + 3exp( 



(B.13) 



CO molecules 

The CO molecule is like the HD molecule and has a permanent dipole moment {^-co = 
0.118 Debye) which is much larger than the dipole moment of the HD molecule (Lang 
1974). The rotational constant is low, Bco ~ 2.78 K, leading to a very low temperature 
for the first transition Ti^ = 5.56 K. 

The expression for the cooling function of CO molecule is obtained in the same way as 
for HD, and we get 
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(B.14) 



and /3q si defined accordingly as in Eqs. (B9) and (BIO). 

H2 molecules 

The H2 molecule differs from HD, since it has no dipole moment and the transitions are 
quadrupolar. We restrict our calculations to the first transition J = -^ J = 2. The 
other transitions can be neglected, since we consider matter temperatures which are much 
lower than the transition temperature of H2 molecules (512 K for the first transition). 
The Einstein coefficients for the spontaneous emission and for the first transition are 
given by (see Kutner 1984): ^^g = 2.44 x 10"" s. 
The relations between the Einstein coefficients are 
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The density of the cosmic background radiation at the transition is: 
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The coUisional coefficients, in the Maxwell-Bolztmann distribution approximation, are 



5^2,0 — ^o,2^^P( ^fr~ ) 
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by analogy with the HD molecule. Because Tj J » T^ we can make the same approxi- 
mations as for the HD molecule, and we find for the populations X^'^ and X^ ^: 
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The probability of radiative de-excitation becomes 
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while the molecular cooling function for the H2 molecule is 
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again /3^2 is defined accordingly as in Eqs. (B9) and (BIO). 



(B.15) 



(B.16) 



B.l. Collision rates 

Since H2 is the most abundant species, the coUisional coefficients are given by 

Cifl - C^S = C^fl =<<J.v> nH2 , (B.17) 

where a ^ lA^ is the collisional cross section, 12^2 the number density of H2 and 
"— ^, with uiH being the hydrogen mass. 
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